Atomistic simulations of adiabatic coherent electron transport in triple donor systems 
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A solid-state analogue of Stimulated Raman Adiabatic Passage can be implemented in a triple 
well solid-state system to coherently transport an electron across the wells with exponentially sup- 
pressed occupation in the central well at any point of time. Termed coherent tunneling adiabatic 
passage (CTAP), this method provides a robust way to transfer quantum information encoded in the 
electronic spin across a chain of quantum dots or donors. Using large scale atomistic tight-binding 
simulations involving over 3.5 million atoms, we verify the existence of a CTAP pathway in a realistic 
solid-state system: gated triple donors in silicon. Realistic gate profiles from commercial tools were 
combined with tight-binding methods to simulate gate control of the donor to donor tunnel barriers 
in the presence of cross-talk. As CTAP is an adiabatic protocol, it can be analyzed by solving the 
time independent problem at various stages of the pulse - justifying the use of time-independent 
tight-binding methods to this problem. Our results show that a three donor CTAP transfer, with 
inter-donor spacing of 15 nm can occur on timescales greater than 23 ps, well within experimentally 
accessible regimes. The method not only provides a tool to guide future CTAP experiments, but 
also illuminates the possibility of system engineering to enhance control and transfer times. 

PACS numbers: 05.60.Gg, 73.63.Kv, 73.21.La 
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I. INTRODUCTION 

The concept of encoding information in the super- 
position of quantum states offers revolutionary ways of 
performing computation and enormous improvement in 
speed and computing power for certain classes of algo- 
rithms Solid-state based quantum computer archi- 
tectures have been the subject of much research due to 
their promise of scalability. Silicon systems are of partic- 
ular interest because of the vast experience of the semi- 
conductor industry in Si electronics and also because Si 
offers a relatively low noise environment for manipulating 
spins. There have been proposals to encode information 
in the nuclear 0] or electronic spin Q of a phosphorus 
donor in Si, the orbital states of a singly ionized two- 
donor molecule Q , the valley-split states of a Si quantum 
well or dot and in gate-confined 2DEGs 

A potentially scalable quantum computer has to in- 
volve complex circuitry of qubits to perform multiple 
levels of error correction and fault-tolerance. In terms 
of the fault-tolerant threshold and defect tolerance such 
architectures will benefit from having separate zones for 
computation and measurement, and hence a mechanism 
for qubit transport 0. Hence, the qubit state typically 
encoded in spin needs nonlocal transport while preserv- 
ing the coherent superposition of the state vectors in the 
Hilbert space - the quantum mechanical equivalent of lo- 
cal bit transfers in traditional computers. 

There already have been several proposals for non-local 



transport of encoded information in solid state quantum 
computers. Skinner et al. [f| proposed a scheme in which 
electrons at an interface between Si and Si02 can be 
shuttled laterally along the surface by appropriate volt- 
age pulses applied to a series of gates. This approach, 
however, requires a high gate density, and is suscepti- 
ble to charge noise and spin-orbit interaction. Other 
approaches have used a chain of coupled harmonic os- 
cillators § or interacting spins [l(J 11 1 in the form of a 
quantum bus. 

In Ref [l2| , a method was proposed to coherently trans- 
port quantum information encoded in the spin of an elec- 
tron across a chain of quantum dots or ionized donors by 
modulating the tunnel barriers between them with volt- 
age pulses. This technique is the solid-state analog of 
Stimulated Raman Adiabatic Passage (STIRAP) from 
quantum optics, and presents a robust population trans- 
fer mechanism by only electrical control. In this process, 
the electron is directly transferred from one end of the 
chain to the other with exponentially suppressed occu- 
pation in the middle of the chain at any point of time. 
This is made possible by adiabatically following certain 
pathways in the eigenspace connecting end states of the 
chain. This technique has been termed coherent tun- 
neling adiabatic passage (CTAP) In addition to 
being a coherent transfer mechanism, this scheme does 
not change the energy of the electron, and is thus ide- 
ally a dissipationless technique. Moreover, this method 
requires gating only donors at the two ends of the chain, 
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and hence can reduce the gate density of the architecture, 
although the three donor CTAP involves no advantages 
in gate density reduction. CTAP was incorporated in the 
bi- linear Si:P architecture design of Ref 0]. 

Previous theoretical works on CTAP have inve stig ated 
its feasibility for transporting sing le atoms [H 01 and 
Bose-Einstein condensates 15|, LL6| ■ CTAP has also been 
recently witnessed using photons in tri ple core optical 
waveguides [13, 18, [Hj]. Recent papers 2(], 21 1 made a 
thorough comparison between the quantum optics and 
solid-state versions of CTAP, highlighting the important 
differences between the two frameworks. Another recent 
work [22| showed that the time of electron transfer in 
CTAP scales as the square root of the number of dots 
in the chain, and the prospects for ion-implanted CTAP 
devices are discussed in Ref [23 1 . 

Here we demonstrate through numerical modeling the 
existence of CTAP in a realistic solid state system taking 
into account the atomistic nature of the underlying semi- 
conductor. Our implicit goal is also to provide incentive 
and guidance to potential CTAP experiments in donors 
or quantum dots. However, we also provide results that 
show the utility of atomistic approaches to the solution 
of demanding time-domain quantum coherent problems. 

The test case investigated here involves a lattice of 28 Si 
atoms with three ionized donors and one bound electron 
under multiple gates, as shown in the schematic of Fig 
1. This is a protoype case for CTAP in a long chain of 
donors as it involves most of the essential physics of the 
many-donor chain. 

CTAP in the quantum optics framework relies on an 
ideal localization assumption, and is well described by a 
three state system 2(J 21 1 . In a solid state system under 
gate bias, this ideal localization assumption is no longer 
valid as the relevant states can admix with higher lying 
excited states. Since we diagonalize the full Hamiltonian 
of the system in an atomistic basis set, the excited states 
are explicitly included in our calculations, enabling us to 
verify whether CTAP can indeed occur in realistic solid- 
state systems once the ideal localization assumption is 
relaxed. 

Large scale atomistic tight-binding device simulations 
perfomed here also enable us to incorporate the Si host 
atoms in the model beyond the effective mass approx- 
imation and hence to include effects due to the full- 
bandstructure of the host. Furthermore, we utilize P 
donor models with valley-orbit interactions and core- 
corrections [24j|, and also use gate potentials obtained 
from commercial Poisson solvers to describe realistic de- 
vices. We are thus able to simulate gate control of tun- 
neling barriers between the donors in the presence of gate 
crosstalk. The simulations carried out in this work are 
some of the most intensive single atom level quantum 
control simulations performed in a realistic solid-state 
system. 

Overall, the technique used here not only gives us 



a powerful simulation tool to model and guide future 
CTAP experiments, but also to show the existence of an 
adiabatic path for solid-state CTAP. This work also sets 
the stage for possible future investigations of the sensi- 
tivity of the adiabatic pathway to donor positioning and 
also of the scalability of the system with increasing num- 
ber of donors. In Section II, we describe a typical CTAP 
device used in this work. Section III explains the con- 
cept of CTAP using a toy model. Section IV outlines the 
tight-binding method and relevant details. In Section V, 
we discuss the three-donor molecular spectrum and tun- 
neling rates. Section VI describes pulse engineering to 
realize CTAP and reports adiabatic transfer times. 
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FIG. 1: Schematic of the device for investigating CTAP. The 
barrier (B) gates modulate the tunneling barriers between the 
donors, while the symmetry (S) gates detune the left and right 
donor eigen states. The simulation has a 3D domain of 60.8 
nm x 30.4 nm x 30.4 nm, and contains about 3.5 million 
atoms. The three donors are labeled L (left), M (middle), 
and R (right), a) Top view b) 3D view. 



II. DEVICE STRUCTURE 

A typical CTAP device design used in this work is 
shown in Fig. 1. The quantum mechanical eigenstates 
are computed in a cubic domain of 60.8 nm x 30.4 nm 
x 30.4 nm with comprising of about 3.5 million Si atoms 
with one electron bound across three ionized P donors. 
The modeling of this device uses an atomistic grid of 3D 
geometry. 

Two symmetry gates and two barrier gates are placed 
above a 5 nm thick oxide layer on top of the Si lattice. 
The function of the barrier (B) gates is to modulate the 
tunnel barriers between the donors. The symmetry (S) 
gates are used to detune the energies of the end donors, 



adding an extra degree of control for the CTAP path- 
way. The gates considered here are of 10 nm width. Due 
to the relative closeness of the gates in this test device, 
there is significant cross-talk between them, making con- 
trol relatively hard for this small device. Gate cross-talk 
in this context means that a typical gate can affect parts 
of the device in addition to its intended function. For 
example, the barrier gate Bi is intended to control the 
tunnel barrier between the left (L) and the middle (M) 
donor only, but in reality it also affects both the barrier 
between the middle and the right (R) donor and the de- 
tuned energies of the end donors, L and R, relative to 
M. In a larger device, the gates will be farther apart, 
reducing the cross-talk effects and easing controllability. 

The donors are buried 15 nm below the oxide, and are 
also placed 15 nm apart from each other in the [100] direc- 
tion. Relatively small device dimensions were chosen as a 
large number of tight-binding simulations had to be per- 
formed for different voltage configurations to zone in on 
the adiabatic path. However, the dimensions used here 
are sufficient to capture all the essential physics of the 
donor chain. Increasing the device domains further only 
leads to an increase in compute times without adding 
new effects to the analysis. The choice of 15 nm distance 
between a donor and an interface is due to the fact that 
a donor located more than 15 nm away from any inter- 
face in this model can be considered as an isolated donor 
in bulk Si free from interface effects (25|. In contrast, 
an ideal CTAP device may have the impurities buried 
deeper and farther apart, and would involve a chain of 
donors in place of the middle donor. 

III. EFFECTIVE 3x3 MODEL OF CTAP 

The concept of CTAP is best described by a simple 
effective model. Assuming three different donor sites, 
and a wavefunction localized in each donor, we can use a 
3x3 Hamiltonian describing the system in this 3-state 
basis. The Hamiltonian H is of the form, 



El t'LM 

H = I thM* Em tMR 
tLR* t M R* E R 



(1) 



where Ei is the on-site energy of i-th impurity, and Uj is 
the tunneling matrix element from impurity i to impurity 
j. We can further simplify the system by assuming the 
ground state of the donors arc aligned in energy so that 
El = Em = En, and arbitrarily set the eigenvalues to 
0. We can also assume only nearest donor coupling by 
setting tLR = 0. The reduced matrix is of the form, 



t LM 
% = I t L M* t M R 
tMR* 



(2) 



The eigenvalues and eigenvectors of this matrix are of 
the form, 
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We initialize the system by localizing the electron 
around the left donor (setting tLM — to prevent hy- 
bridization). Looking at Eqs. 6-8, we see that this con- 
figuration corresponds to the second dressed state | \I> 2 ) - 
In a real system, setting tLM = involves raising the 
tunneling barrier between L and M by applying a more 
negative bias to Bi, while adjusting other gates to com- 
pensate for cross-talk. In practice, this also means that 
\tLu\ <C \tMB.\ so that the barrier between M and R is 
lower than that between L and M. At some point later we 
set the tunnel barriers tLM = tMR, at which point \H>2) 
represents a superposition of the L and R localized states 
(with no appreciable population in state M). Finally, we 
reduce tMR smoothly to zero leaving the system localized 
in state R. Following such a sequence allows us to adia- 
batically evolve the system from a L localized state to a 
R localized state, without populating the central donor. 
This adiabatic pathway is described in the ideal limit 
by simply following state ^2) from initial to final state. 
During the entire transfer process, E 2 is held fixed at 0, 
resulting in no change of electron energy and no acquisi- 
tion of any dynamical phase. Since this is an adiabatic 
transfer process at a very low temperature, the electron 
always occupies the state in which it started. 

The pulsing sequence described above is quite counter- 
intuitive in nature. If we are trying to transfer the elec- 
tron from L to R, then an intuitive pulsing sequence 
would involve lowering the barrier between L and M first, 
and then the barrier between M and R. This will trans- 
port the electron first to the middle well, and then to the 
right well, very much like a bucket-brigade device. How- 
ever, if the barriers are modulated in the reverse order 
such that the barrier between M and R is lowered first 
and then that between L and M according to the CTAP 
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protocol, then the electron is transferred directly from L 
to R in a much more robust fashion (l2| in terms of pulse 
control over tunnel rates. A signature of the CTAP pro- 
tocol is an exponentially suppressed occupation at the 
middle donor. 



IV. METHOD 

The tight binding model employed in this work is the 
20 band sp 3 d 5 s* spin model with nearest-neighbour in- 
teractions. This model incorporates spin inherently in 
the basis by duplicating the 10 spatial orbitals per atom 
for up and down spins. Spin-orbit interactions of the 
host are also included by onsite p-orbital spin-orbit cor- 
rections [26| . The model parameters were optimized by a 
genetic algorithm with appropriate constraints to repro- 
duce the important features of the bulk bandstructure 



of the host [23, 12a I2S1- The p donors were modeled by 
Coulomb potentials screened by the dielectric constant 
of Si. At the donor site, a cut-off potential Uq, was used, 
and its value optimized so that the ground state binding 
energy of -45.6 meV was obtained for a donor in bulk Si. 
In this model, the magnitude of Uo reflects the strength 
of the valley-orbit interaction responsible for lifting the 
six-fold degeneracy of the Is manifold of the impurity. 
It was shown in an earlier work [24[ that the splitting 
between the singlet, triplet and doublet components of 
the Is manifold increase with the magnitude of Uo- This 
semi-emprical highly optimized technique was able to re- 
produce the full single donor spectrum very accurately. 

The electrostatic gate potentials were obtained from 
a commercial Poisson solver, and was then interpolated 
into the atomistic grid for the tight-binding simulations. 
Due to the large number of computer-intensive simula- 
tions required to home in on the adiabatic path, it was 
not possible to generate each time the total potential 
profile of all the gates taken together. Instead, we gen- 
erated the potential profiles of each gate separately, and 
assumed the net potential can be obtained from the su- 
perposition principle. Some non-linear behavior is ex- 
pected to arise from the additional fringing fields near 
the gates when multiple gates are turned on at the same 
time. However, such non-linear behavior is expected only 
to add small voltage corrections. Furthermore, the ba- 
sic principle of tuning the tunnel barriers between the 
donors by barrier gates to realize the adiabatic pathway 
remains unchanged. 

Closed boundary conditions with a model of dangling 
bond passivation was used to model the interfaces. The 
full Hamiltonian of about 3.5 million atoms including the 
four gate potentials was solved by parallel Lanczos and 
Block Lanczos algorithm to capture the relevant eigen- 
values and wave functions. Typical computation time 
for 6 states was 7 hours on 40 processors [30(. Although 
CTAP is a time dependent problem requiring transient 



voltage pulses, it can be analyzed by snapshots of the 
wave- functions at different biases obtained from the time 
independent Schrodinger equation. 

The tight-binding method under the hood of the Nano- 
electronic Modeling Tool (NEMO-3D) [H, was used 
previously to verify the Stark shift coefficients of the hy- 
perfine interaction of the donor spin [3~J ] with respect to 
ESR measurements [32|. The method was also used to 
compute orbital Stark shifts of an As donor in Si close to 
the oxide barrier, and could explain energy level measure- 
ments from transport experiments in commercial FIN- 
FETs [33j . The same method was successfully applied to 



investigate valley-splitting with alloy-disorder and step 
roughness in Si quantum wells [23|, and was also used 
to model quantum dots for optical communication wave- 
lengths [35j. 
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FIG. 2: The energies of the 4 lowest single electron states of 
the 3P 2+ donor molecule as a function of donor separation 
in [100]. The energies are shown relative to the conduction 
band minima. The energy differences between the states cor- 
respond to tunnel couplings between the donors. 



V. THREE-DONOR MOLECULE AT ZERO GATE 

BIAS 



A single phosphorus donor in bulk Si has a ground 
state (GS) of Ai symmetry at -45.6 meV relative to the 
conduction band (CB) minima. Above this, there is an 
orbital triplet manifold of T2 symmetry at -33.9 meV, 
and an orbital doublet manifold of Ei symmetry at -32.6 
meV [37J. When three ionized donors are located close- 



by, coupling between the wells produce molecular states 
that may span over the whole chain. Fig. 2 shows the 
energies of the four lowest states of the donor molecule 
(3P 2+ ) at zero gate bias as a function of separation dis- 
tance along [100]. The energy differences between these 
states are proportional to the tunnel barriers. The sep- 
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aration distance between the donors is incremented in 
equal steps so that the L and R are always equidistant 
from M. 

It has been shown that the tunnel coupling for a two- 
donor (2P + ) charge qubit exhibits oscillatory behaviour 
with relative donor separations along [110] and [111] [38| . 
The exchange coupling between the donor electrons of a 
two-donor molecule (2P) has also been shown to exhibit 
oscillations with donor positions, a consequence of phase 
pinning of the Bloch functions at the donor sites [39I ]. 
This has posed some controllability issues for two-qubit 
operations of the Si:P based qubits 

USUI, although 

individual qubit characterization is likely to resolve the 
problem. It would be interesting to investigate if the 
CTAP protocol is also susceptible to radial and angular 
donor misalignments, but this goes well beyond the scope 
of our present work. 
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FIG. 3: Particular potential landscapes (a4 and b4) and the 
corresponding molecular states (al-3 and bl-3) of the triple- 
donor device. The left column shows the donor molecule at 
zero gate bias. The right column shows the donor molecule 
at the middle of the adiabatic path. Note the almost zero 
electron density at the middle donor in the 2nd eigenstate. 

In a triple Coulomb well generated from the superpo- 
sition of three isolated Coulomb potential, the middle 
well is deeper than the left and right wells. As a result, 
the ground state of the system will have electron density 
at the middle donor. The L and R wells are essentially 
at the same energy, and form nearly degenerate bonding 
and anti-bonding states. Top row of Fig 3 shows the to- 
tal device potential including the donors and the gates, 



while the remaining rows show the three lowest wave- 
functions of 3P 2+ . The left column of Fig. 3 portrays 
the donor molecule at zero gate bias. The wavefunctions 
conform to the symmetries described above. At low gate 
bias, these lowest three states arise from linear combina- 
tions the single Ai states from each donor. Similarly, the 
closely spaced orbital triplet and doublet manifolds from 
each impurity will interact to form 15 excited molecular 
states. Since there is a gap of about 11 meV between the 
orbital singlet and the orbital triplet manifold of a single 
impurity, there will also be a modest energy gap between 
the lowest three molecular states and the higher states. 
However, this gap is likely to decrease as donor separa- 
tions decrease, and we need to make sure that the mani- 
fold of the lowest three states is sufficiently isolated from 
the higher states for ideal CTAP operation. Fig. 2 also 
shows that the centrally occupied state Ei approaches 
E2 and E3 as donor separation increases and the system 
moves towards the isolated donor regime, in which all the 
three Ai states are degenerate. E4 is the lowest molecu- 
lar state arising from the T2 manifold. Even if the gate 
lengths can be scaled down to nanometers, donor sepa- 
rations less than 10 nm are not desirable as the E4 state 
approaches the lowest three manifold. However, donor 
separations of 15 nm or more seem reasonable for CTAP. 



VI. GATE PULSING TO REALIZE CTAP 

To obtain the wavefunction symmetries for CTAP de- 
scribed in Section III, the molecular states needed to be 
aligned close to each other in energy. This was achieved 
by applying negative biases to the barrier gates so that 
the middle donor is pulled up in energy close to E2 and 
E3. This also lowers the effective tunnel barriers be- 
tween the donors causing more hybridization between the 
states. In the right column of Fig 3, we applied -0.11 V 
to each of the barrier gates, and compensated for gate 
crosstalk by applying a small bias to the symmetry gate 
Si. The wavefunctions indeed correspond to the symme- 
tries identified from the effective 3x3 model. However 
there is only a limited range of barrier gate biases where 
this happens. In the simulations, we noticed a barrier 
gate bias window between -0.1 V and -0.12 V where the 
states are strongly interacting to produce the molecular 
symmetries we seek. 

With small barrier gate modulation around -0.11 V, 
we were able to obtain other points on the adiabatic tra- 
jectory. In Fig 4, we show the eigenvalues (4a) and the 
barrier Bi and B2 gate voltages (4b) at 10 points of the 
adiabatic path. We observe anti-crossing between the 
states Ei and E2 at the middle of the path in Fig 4a. The 
energy gap between the lowest two states approaches a 
non-zero minimum of 29 /ieV at the middle of the path 
at point 5. From this minimum gap, we can roughly esti- 
mate the electron transfer timcscalc as r = — = 0.023 
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FIG. 4: (a) Evolution of the lowest three eigenvalues with the 
CTAP gate pulse. The gap between Ei and Ei is minimum 
at point 5. (b) The Bi and B2 gate voltages at the 10 points 
of the adiabatic path. At all points, the Si and S2 gates are 
held fixed at -0.01325 V and V respectively to compensate 
for cross-talk. 



ns for a distance of 30.4 nm. This minimum energy gap 
depends both on barrier gate voltages and donor sepa- 
rations, and can be engineered to achieve faster transfer 
times. 

Fig 5 shows the electron localization corresponding to 
the first excited state at various stages of the electron 
transfer. The population gradually diminishes in the left 
donor and reappears in the right, with minimal leakage 
to the middle donor. This verifies for the first time that 
CTAP can be observed in a realistic solid-state system 
such as an interacting donor chain. Unlike the solution 
of the effective 3x3 model outlined in Section III, the 
more realistic simulations show that the CTAP state has 
a non-zero and time-varying energy. However, this is 
still essentially the same CTAP protocol described by 
the effective 3x3 model. 



VII. CONCLUSION 

We have demonstrated the possibility of CTAP in a 
triple donor chain through precise numerical modeling. 
Originally developed in the quantum optics framework, 
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FIG. 5: The 1st excited state under a voltage sequence that 
transfers the electron from the left to right donor without 
middle occupation. The voltage configurations involve mod- 
ulations of Bi and B2 gate voltages, while keeping Si and S2 
fixed. Only 5 points of the adiabatic path are shown here. For 
the barrier gate voltages at these points, refer to Fig 4(b). 



CTAP relies on an ideal localization assumption and a 
few state based description, which are generally not met 
in solid-state systems. We have shown that a realistic 
solid-state system can still exhibit CTAP when the few 
state description is abandoned and a Schrocdinger wave 
description is used with many molecular states consid- 
ered in the calculations. Despite controllability prob- 
lems due to gate crosstalk in a small device and band- 
structure effects of the host material, it is possible to find 
adiabatic trajectories that define CTAP. The large scale 
highly optimized quantum mechanical device simulations 
done here not only show the existence of adiabatic path- 
ways in a triple donor system, but also helps to devise 
a technique to model and guide potential experiments. 
The results enable us to estimate typical adiabatic trans- 
fer timescale of 0.023 ns for this device with left and right 
impurities separated by 30.4 nm. Since CTAP presents a 
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robust and coherent method to transport electronic spin 
in a quantum circuitry, experimental demonstration of 
CTAP in solid-state systems should be sought after. 

Although the three-donor case serves well as a test, 
the real benefit of CTAP will be evident in a long donor 
chain. Under suitable gate pulses, the donor electron 
can be transported from one end of the chain to another, 
carrying along with it the quantum information encoded 
in its spin. Realizing such a system will indeed help to 
solve some of the critical information transport problems 
in solid state quantum computing architectures. It is 
therefore necessary to investigate scalability of the adi- 
abatic pulsing scheme to increasing number of donors. 
Further studies need to be undertaken to investigate the 
sensitivity of the adiabatic path to relative donor posi- 
tioning and also to investigate spin-orbit coupling effects 
at various stages of the transfer. 
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